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Abstract 

In the minimal supersymmetric standard model (MSSM), the neutralinos, the spin- 
1/2 Majorana superpartners of the neutral gauge and Higgs bosons, are expected to 
be among the light supersymmetric particles that can be produced copiously at future 
high-energy colliders. We analyze two-body neutralino decays into a neutralino plus 
a Z boson or a lightest neutral Higgs boson h, allowing the relevant parameters to 
have complex phases. We show that the two-body tree-level decays of neutralinos 
are kinematically allowed in a large region of the MSSM parameter space and they 
can provide us with a powerful probe of the Majorana nature and CP properties of the 
neutralinos through the Z-boson polarization measured from Z-boson leptonic decays. 
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1 Introduction 



The search for supersymmetry (SUSY) is one of the main goals at present and future colliders 
PQ as SUSY is generally accepted as one of the most promising concepts for physics beyond 
the Standard Model (SM) |2j. A special feature of SUSY theories is the existence of the 
neutralinos, the spin-1/2 Majorana superpartners of the neutral gauge bosons and Higgs 
bosons. In the MSSM, the neutralinos are expected to be among the light supersymmetric 
particles that can be produced copiously at future high-energy colliders Once several 
neutralino candidates are observed at such high-energy colliders, it will be crucial to establish 
the Majorana nature and CP properties of the neutralinos. In this light, many extensive 
studies of the general characteristics of the neutralinos in their production and decays 
ElEj as well as in the selectron pair production [E] at e + e~ and/or e~e~ linear colliders have 
been performed. 

In the present work, we analyze two-body tree-level decays of neutralinos into a neu- 
tralino plus a Z boson or a lightest neutral Higgs boson h in order to probe the Majorana 
nature of the neutralinos and CP violation in the neutralino system. A comprehensive analy- 
sis of the two-body decays of neutralinos as well as charginos was given previously in Ref . jH] • 
We however note that a rather light Higgs boson mass was assumed and no Z boson po- 
larization was considered in the previous work. One powerful diagnostic tool in the present 
analysis is Z polarization, which can be reconstructed with great precision through Z-boson 
leptonic decays, Z — > in particular, with I = e,fi. 

It is possible that due to the masses of the relevant particles, no two-body tree-level 
decays are allowed, in which case the dominant decays would consist of three-body tree- 
level PO] or two-body one-loop decays [TT]. However, a sufficiently heavy neutralino can 
decay via tree-level two-body channels containing a Z or h with its mass less than 135 
GeV in the context of the MSSM If some sfermions are sufficiently light, two-body 

tree-level decays of neutralinos into a fermion and a sfermion may be also be important. 
However, neutralinos heavier than the squarks will be extremely difficult to isolate at hadron 
colliders, because the squarks and gluinos are strongly produced and they decay subsequently 
into lighter neutralinos and charginos. On the other hand, at e + e~ colliders, squarks and 
sleptons, if they are kinematically accessible, are fairly easy to produce and study directly. 
With these phenomenological aspects in mind, we assume in the present work that all the 
sfermions are heavier than (at least) the second lightest neutralino Xt Then, we investigate 
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the MSSM parameter space for the two-body tree-level decays of the neutralino x® an d show 
how the Majorana nature and CP properties of the neutralinos can be probed through the 
two-body decays x\ ~~ *" Xi Z, once such two-body decays are kinematically allowed. 

The paper is organized as follows. Section |2] is devoted to a brief description of the mix- 
ing for the neutral gauginos and higgsinos in CP-noninvariant theories with non-vanishing 
phases. In Sec. El after explaining the reconstruction of Z-boson polarization through the 
Z decays into two-lepton pairs, we present the formal description of the (polarized) decay 
widths of the two-body neutralino decays into a lightest neutralino x\ pl us a % boson or a 
lightest Higgs boson h with special emphasis on the polarization of the Z boson. In Sec. 01 
we first investigate the region of the MSSM parameter space where the two-body neutralino 
decays are allowed and discuss the dependence of the branching ratios and decay widths on 
the relevant SUSY parameters. Then, we give a simple numerical demonstration of how the 
Majorana nature and CP properties of the neutralinos can be probed through the two-body 
decays X2 ~* Xi Finally, we conclude in Sec. 



2 Neutralino Mixing 



In the MSSM, the mass matrix of the spin- 1/2 partners of the neutral gauge bosons, B and 
W 3 , and of the neutral Higgs bosons, H® and H®, takes the form 



M 



N 



( M 1 


-nizcpsw 
\ mzspsw 





M 2 
nizcpcw 



-m z cps w 
m z cpc w 


-fj, 



m z sps w \ 
-m z spc w 
-H 

) 



(1) 



in the {B, W 3 , H®, H$} basis. Here M x and M 2 are the fundamental SUSY breaking U(l) 
and SU(2) gaugino mass parameters, and fi is the higgsino mass parameter. As a result 
of electroweak symmetry breaking by the vacuum expectation values of the two neutral 
Higgs fields v\ and v 2 {sp = sin/3, cp = cos/? where tan/5 = t^/fi), non-diagonal terms 
proportional to the Z-boson mass mz appear and the gauginos and higgsinos mix to form 



the four neutralino mass eigenstates x 



1-4), ordered according to increasing mass. 



In general the mass parameters Mi, M 2 and \x in the neutralino mass matrix (0) can be 
complex. By re-parameterization of the fields, M 2 can be taken real and positive, while the 
U(l) mass parameter Mi is assigned the phase $1 and the higgsino mass parameter \x the 
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phase For the sake of our latter discussion, it is worthwhile to note that in the limit 
of large tan/? the gaugino-higgsino mixing becomes almost independent of tan/3 and the 
neutralino sector itself becomes independent of the phase $ M in this limit. 

The neutralino mass eigenvalues = m^q (i = 1-4) can be chosen positive by a suitable 
definition of the mixing matrix N, rotating the gauge eigenstate basis {B, W 3 , H®, H®} to 
the mass eigenstate basis of the Majorana fields: N*M. ^ N< = diag(m^o, m^o, m^o, m^,o). In 
general the mixing matrix iV involves 6 non-trivial angles and 9 non-trivial phases, which 
can be classified into three Majorana phases and six Dirac phases [4 . The neutralino sector 
is CP conserving if /i and Mi are real, which is equivalent to vanishing Dirac phases (mod 
7r) and Majorana phases (mod vr/2). Majorana phases of ±7r/2 do not signal CP violation 
but merely indicate different intrinsic CP parities of the neutralino states in CP-invariant 
theories [Zj. 



3 Two— Body Decays of Neutralinos 

Before describing the two-body decays x® ~ * Xj Z in detail, we explain how to reconstruct 
the Z polarization through the lepton angular distributions of the Z-boson leptonic decays, 
Z — ► l~l + , particularly with I = e, \x. In the rest frame of the decaying Z boson, which can 
be reconstructed with great precision by measuring the lepton momenta, the lepton angular 
distributions are given by 

1 dT[Z{±) -> i+r] 3 



T[Z^m-] d cos 6i I 

1 dT[Z(0) -> _ 3 , 



cos $i ± 2 £i cos Q\ 



for the Z-boson helicities, ±1 and 0, respectively, where £j = 2viai/(vf + af) ~ —0.147 with 
v i — s w — 1/4 and a\ = 1/4, and d\ is the polar angle of the l~ momentum with respect to 
the Z boson polarization direction. Here, the decay width T[Z — > is the average of 

three polarized decay widths, 

r[z -> i + r] = ^ { r[z(+) -> i + r] + r[z(o) -> i + r] + r[z(-) -> i+r) } . (3) 

We emphasize that the three polar-angle distributions (j2J can be determined without know- 
ing the full kinematics of the decay xl ~^ Xj Z '■ m contrast, the distributions involving 
the interference of the amplitudes with different Z helicities are always accompanied with 
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azimuthal angle dependent terms. As the lightest neutralino Xi assumed to be the light- 
est SUSY particle (LSP) always escapes detection, the kinematics of the two-body decay 
X.1 ~" * Xj Z cannot be fully reconstructed so that the azimuthal-angle dependent distribu- 
tions are not fully available. 

The decay width of the decay xi 
reads 



Xj Z producing a Z boson with its helicity, ±1 or 0, 



xi 



X.i 
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(4) 



16irmf 

respectively, where the asymmetries A^ and A? are defined in terms of the vector and 
axial-vector couplings V and A of the Z boson to the neutralino current as 



A 



N 



\A\ 



\V\ 2 + \A\ 



A 
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2Me{VA* 
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(5) 



771; 



m 3) 



rrir 



with gz = gj cos8w and the kinematical factor Xz = [(rrii + rrij) 2 — m 2 z }[( 
Combining the leptonic Z-boson decay distributions (0) with the polarized decay widths 
we obtain the correlated polar-angle distribution: 
dT 



d cos 9i 



.B[z-.zi]{(r[x?-,jgz(+)] + r[x»- 

+ 2{v[x1^x]Z{+)]-V[x\ 
+ 2r[x?-$Z(0)] sin 2 ^}. 



xjZi-)}) (1 + cos 2 ^ 
-X? £(-)]) a cosfl. 



(6) 



Consequently, each polarized decay width can be extracted from the correlated polar-angle 
distribution by projecting out the distribution with a proper lepton-polar angle distribution. 

The explicit forms of the vector and axial-vector couplings V and A in Eq. (jSJ) are given 
in terms of the 4x4 neutralino diagonalization matrix N in the MSSM by 



V = --Qm (N j3 N* 3 - N jA N& , A = (n j3 n* 3 - n^n*,) . (7) 

Note that the vector coupling V is pure imaginary and the axial-vector coupling A is pure 
real. This characteristic property of the Z-x1-X°j coupling due to the Majorana nature of 
neutralinos leads to one important relation between the polarized decay widths with the 
Z-boson helicities, ±1: 



r[x?-xSz(+)] = ml^x*z{-)\ 



(8) 
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which is valid even in the CP non-invariant theory. This relation can be checked by measur- 
ing the forward-backward polar-angle asymmetry of the correlated polar-angle distribution 
(JOJ). However, because of the small analyzing power ~ —0.147, it will be necessary to have 
sufficient large number of decay events to measure the asymmetry with good precision. In 
addition to the relation (jSJ), the relative intrinsic CP parity of two neutralinos in the CP 
invariant theory can be determined by measuring the ratio of the longitudinal decay width 
to the transverse decay width, which satisfies 

LT "r[x°-x^(+)] + r[x°-x^(-)] m% ' w 

for the even/odd relative intrinsic CP parity with V = 0/A = 0, i.e. An = ^Ij respectively. 
In the CP non-invariant theory, both the vector and axial couplings are in general non- 
vanishing, leading to the value of the asymmetry different from ±1. Therefore, any 
precise measurements of the asymmetry A^ will provide us with an important probe of 
CP violation in the neutralino system under the assumption that the neutralino masses are 
measured with good precision, independently of the decay modes. 



Next, we give the decay formulas into final states containing a lightest neutral Higgs 

o 
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boson h. The explicit form of the decay width of the decay Xi ~ > X? h is written as 



Xi - *S h \ = £rK \\ S \ 2 ( K + m,) 2 - ml) + \P\ 2 (K " m,) 2 - ml)} , (10) 



with the kinematical factor = [(m.; + mj) 2 — m^] [(rrij — mj) 2 — m^]. The scalar and 
pseudoscalar couplings, S and P, of the Higgs boson h to the neutralino current are defined 
in terms of the mixing matrix iV as 

S = ^ 9£e \(N j2 - twNjiXsaNx + c a N u ) + (i <-> j)) , 

P = [(N j2 - twN^SaNa + c a N iA ) + (i <-> j)} , (11) 

where tw = tan.9w, c a = cos a and s a = sin a for the neutral Higgs mixing angle a. If the 
charged Higgs boson mass in the MSSM is very large, 

c a — > sin ft, s a ^—cosP, (12) 

This decoupling approximation of the cosine and sine of the mixing angle a is very good if 
the charged Higgs mass is larger than twice the Z boson mass ^Sj- For the sake of discussion, 
we take the decoupling limit in the present work. 
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4 Numerical Analysis of Two— Body Decays 



In some SUSY scenarios, the lightest neutralino Xi is the LSP and the second lightest neu- 
tralino x\ among the other three neutralino states are expected to be lighter than sfermions 
and gluino [3]. Then, the two-body decays x\ ~ y Xi Z or h as we U as the two-body decays 
of the heavier neutralinos x® 4 [9 wm constitute the major decay modes of the neutralinos, 
respectively, once the two-body tree-level decay modes are kinematically allowed. In the fol- 
lowing numerical analysis we will ignore all other modes except for the two-body tree-level 
decays of the neutralinos. 

4.1 Branching ratios 

For the branching ratio calculations for the two-body decays x[> — > Xi Z/h, we assume that 
all the SUSY parameters are real, M\ is related to M 2 by the gaugino mass unification 
condition \M\\ = (5/3) tan 2 9yyM 2 ~ 0.5 M 2 and the Higgs boson mass rrih is 115 GeV. In 
addition, we assume that the MSSM Higgs system is in the decoupling regime so that the 
characteristics of the lightest Higgs boson h is similar to the SM Higgs boson to a good 
approximation [T3j . 

Figure ^ shows the regions of the decays of x 2 011 the {/i, M 2 } plane. In the region 
denoted by "three-body decays", no two-body modes are kinematically allowed. In the 
a Z region" (red-colored), only the two-body decay into a Z is allowed and in the u Z/h 
region", both the two-body decays X2 ~^ Xi^/h are allowed. We divide the u Z/h region" 
into three parts, according to B\Z\ < 10%, 10% < B\Z] < 20% and B[Z] > 20%. (Here, 
B[Z] = Br[^2 ~~ > Xi Z]-) m addition, as a reference, the region excluded by the experimental 
bound on the lighter chargino mass m-± > 104 GeV is displayed by the blue-hatched 
region. 

We first note that, if M 2 < 2m z, the mass difference m^o — m^o is less than mz for all 
fj, and the mass difference is very small for M\, M 2 . So, as clearly shown in Fig. Q the 

two-body decay x° ~^ Xi Z is allowed only when 2m z ^ M 2 < 2|//| under the assumption 
of the gaugino mass unification condition. In addition, we find from the figure that for the 
two-body decays the magnitude of \x is required to be larger than about 270 GeV and that, 
once the two-body Higgs mode x\ ~* Xxh 1S open kinematically, this two-body decay mode 
dominates in most of the Z/h region. The region where the decay x 2 Xi Z is appreciable 
is not symmetric between positive and negative /x in the Z/h region. The branching ratio 
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Figure 1: Three distinct regions of x® decay are exhibited as a function of M 2 and /i, 
assuming that M 2 and \x are real. In the region denoted by "three-body decays", no two- 
body modes for the neutralino x 2 except for the loop-induced two-body radiative decays are 
kinematically allowed. In the "Z region" (red-colored), only the two-body decay into a Z is 
allowed and in the u Z/h region", both the two-body decays are allowed. The u Z/h region" is 
divided into three parts, according to B[Z] < 10%, 10% < B[Z] < 20% and B[Z] > 20%. For 
reference, the exclusion region by the experimental bound on the lighter chargino mass bound 
m x ± — 104 GeV is displayed by the blue-hatched region. In this numerical illustration, we 
set tan f3 = 10. 



B[Z] is significant only in a small area of the positive \x region, but in a large area of the 
negative /i region. 

On the other hand, we find numerically that, for the heavier neutralinos \\ 4 , the {M 2 , \x\ 
region for the two-body decays x*] 4 — ► Xi Z/h expands drastically. A large region with small 
\n\ but large M 2 as well as with small M 2 but large \fi\ also allows for the two-body decays of 
the heavier neutralinos, X34 ~~ Xi Z/h. Only in the wedge-shaped band region of the width 
of about 100 GeV around the line satisfying the relation M 2 ~ 2\fi\ no two-body decays for 
the heavier neutralino X3 are allowed, while the heaviest neutralino can still decay into Xi 
and Z/h in the (almost) entire parameter space, possibly except for the region excluded by 
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the experimental lighter chargino mass bound. 

Consequently, for most of the parameter space of the MSSM the decays of the two heavier 
neutralinos are dominated by two-body tree-level processes of which the final state consists 
of a Z boson or h boson together with one of the lighter neutralinos, or a W boson and one 
of the charginos. Furthermore, the two-body decays of the second lightest neutralino x 2 can 
be significant in a large region of the parameter space of the MSSM. 

M 2 =250 GeV, tan(3=10 

5 i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i 




1000 



\l [GeV] 

Figure 2: The dependence of the decay width T[x 2 Xi Z] on the higgsino mass parameter 
fi, assuming that /i is real and Mi = (5/3)t^/M 2 . For this numerical illustration, we set 
M 2 = 250 GeV and tan/? = 10. 

In addition to the branching ratios, it is also crucial to analyze the absolute size of the 
decay width r[x 2 — * Xi Depending on the values of the relevant couplings, the two-body 
decay widths could be smaller than the three-body decay widths involving virtual sfermion 
exchanges, unless the sfermions are too heavy. We exhibit in Fig. El the dependence of the 
decay width T[\ 2 ~ > Xi^] on the higgsino mass parameter //, assuming again that \x is 
real and taking M x = (5/3) tan 2 9 W M 2 , M 2 = 250 GeV and tan/3 = 10. The decay width 
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decreases rapidly with increasing This is because the couplings of the Z boson to the 
neutralino current are governed by the higgsino components of the neutralinos (see Eq. (JJJ)) 
so that the Z-j^-X^ couplings are strongly suppressed for large Therefore, for large 
some three-body decays could be more dominant than the two-body decays. 

1 
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Figure 3: The dependence of the ratio T CP (red solid line) and the asymmetry An (blue 
dot-dashed line) on the CP phase $i for the set of real parameters, {tan/5 = 10, M 2 = 
250 GeV, = 500 GeV} . The phase $ M is set zero in this numerical illustration figure, 
as the ratio T CP and the asymmetry An are found to be insensitive to the phase $ M for the 
given specific set of real parameters. 

4.2 A probe of CP violation 

In the previous subsection, we restrict ourselves to the CP invariant case with real param- 
eters, as the qualitative results obtained from the CP-even quantities are not expected to 
change so significantly even if the parameters are complex. But, the parameters Mi and 
H are in general complex so that it is important to check whether they indeed have com- 
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plex phases or not. The existence of the complex phases in the neutralino system, which in 
general cause CP violation, can be established by the measurements of the ratio 

Tl LT - (mf + m 2 )/m| 

2m imj /m 2 z ' 1 } 

with TZlt defined in Eq. Q as well as the neutralino masses rrii and rrij. The ratio T CP is 
— 1 (+1) in the CP invariant theory for the positive (negative) relative intrinsic CP parity 
of the neutralinos, xl an d Xji taking part in the decay Xi ~~ * Xj respectively 

To explicitly show the dependence of the ratio T CP on the CP phases $i and we 
chose a specific set of real parameters {tan/? = 10, M 2 = 250 GeV, = 500 GeV} as a 
simple numerical example with \M\\ = (5/3) tan 2 ^M 2 , while varying the phases $x an d 
$ M . Numerically, we find that the ratio T CP is insensitive to the phase $ M . This is because 
the phase dependence is always accompanied with sin 2/3 = 2tan/5/(l + tan 2 /3) « 0.2 for 
tan (3 = 10, which is already small, and the higgsino components of the neutralinos are small 
for large So, we show in Fig. |3]the dependence of the ratio T CP (red solid) as well as the 
asymmetry An (blue dot-dashed) only on the phase $i for the real parameter set for one 
fixed value of $ M = 0. Clearly, in the CP-invariant case with $! = 0, tt or 2tt, the absolute 
magnitude of the ratio T CP as well as the asymmetry An is 1, but it is different from 1 in the 
CP non-invariant case. In the given real parameter set, we find that the ratio T CP is quite 
sensitive to the phase $i near $i = it, while it is not so sensitive to the phase near $x = 
and 27T. 

In the present work the analysis for probing the Majorana nature and CP violation in 
the neutralino system has been carried out at the tree level. However, it will be important 
to include loop corrections to the two-body tree-level decays because, if they are small, 
the tree-level CP violation effects might be diluted by loop-induced CP violation effects 
originating from other sectors of the MSSM. 



5 Conclusions 

For a large portion of the MSSM parameter 
X2 as we U as the heavier neutralinos X34 
which the final state consists of a Z or a 
neutralino x?) assumed to be the lightest 



space, the decay of the second lightest neutralino 
could be dominated by two-body processes in 
lightest Higgs boson h together with a lightest 
supersymmetric particle. The main conclusion 
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of the present work is that, unless the two-body decay Xi ~ > Xj Z is strongly suppressed, 
the Z polarization, which can be reconstructed through great precision via the leptonic Z— 
boson decays Z — > provides us with a powerful probe of the Majorana nature of the 

neutralinos and CP violation in the neutralino system. 



Acknowledgments 

The work of SYC was supported in part by the Korea Research Foundation Grant (KRF- 
2002-070-C00022) and in part by KOSEF through CHEP at Kyungpook National University 
and the work of YGK was supported by the Korean Federation of Science and Technology 
Societies through the Brain Pool program. 



References 

[1] TESLA Technical Design Report, Part: III Physics at an e + e~ Linear Collider, eds. 
R.-D. Heuer, D. Miller, F. Richard and P. Zerwas, DESY 2001-011 |hep-ph/0106315] ; 
T. Abe et al. [American Linear Collider Working Group Collaboration], "Linear 
collider physics resource book for Snowmass 2001. 2: Higgs and supersymmetry 
studies", in Proc. of the APS/DPF/DPB Summer Study on the Future of Parti- 



cle Physics (Snowmass 2001) , ed. N. Graf, hep-ex/0106056; K. Abe et al, JLC 



Roadmap Report, presented at the ACFA LC Symposium, Tsukuba, Japan 2003, 



http://lcdev.kek.jp/RMdraft/; G. Guignard (ed.), A 3-TeV e + e linear collider based 



on CLIC technology, CERN-2000-008. 

[2] H.P Nilles, Phys. Rept. 110 (1984) 110; H. Haber and G. Kane, Phys. Rept. 117 (1985) 
75. 



[3] See, for instance, B.C. Allanach et al., Eur. Phys. J. C 25 (2002) 113 |hep-ph /0202233 . 



[4] SY. Choi, J. Kalinowski, G. Moortgat-Pick and P.M. Zerwas, Eur. Phys. J. C 22 563 
(2001) lhep-ph/0108"TT7| , ibid. C 23 769 (2002) hep-ph/ 0202039"! ; J. Kalinowski, Acta 
Phys. Polon. B34 (2003) 3441 l hep-ph/0306272l ; SY. Choi, |hep-ph/0308060| 

[5] S.T. Petcov, Phys. Lett. B139 (1984) 421; S.M. Bilenky, E.Kh. Khristova and N.P 
Nedelcheva. Bulg. J. Phys. 13 (1986) 283; G. Moortgat-Pick and H. Fraas, Phys. Rev. 



12 



D 59 (1999) 015016 |hep-ph/9708481|; G. Moortgat-Pick, H. Fraas, A. Bartland W. Ma- 



jerotto, Eur. Phys. J. C 9 (1999) 521 [Erratum-ibid. C 9 (1999) 549] |hep-ph/9903220 ; 
G. Moortgat-Pick, A. Bartl, H. Fraas and W. Majerotto, Eur. Phys. J. C 18 (2000) 379 
|hep-ph/0007222| ; M.M. Nojiri, D. Toya and T. Kobayashi, Phys. Rev. D 62 (2000) 
075009 |hep-ph/0001267| . 

[6] S.M. Bilenky, N.D. Nedecheva and S.T. Petcov, Nucl. Phys. B247 (1984) 61; S.T. Pet- 
cov, Phys. Lett. B178 (1986) 57; N. Oshimo, Z. Phys. C21 (1988) 129; Y. Kizukuri and 
N. Oshimo, Phys. Lett. B249 (1990) 449; S.Y. Choi, H.S. Song and W.Y. Song, Phys. 
Rev. D 61 (2000) 075004 | hep-ph/9907474| ; V.D. Barger, T. Falk, T. Han, J. Jiang, T. 
Li and T. Plehn, Phys. Rev. D 64 (2001) 056007 |hep-ph/010TT06| ; A. Bartl, H. Fraas, 



O. Kittel and W. Majerotto, hep-ph/0308141 A. Bartl, T. Kernreiter and O. Kittel, 



|hep-ph/ 0309340; S.Y. Choi, M. Drees, B. Gaissmaier and J. Song, |hep-ph/0310284| 

[7] J. Ellis, J.M. Frere, J.S. Hagelin, G.L. Kane and S.T. Petcov, Phys. Lett. B 132 (1983) 
436; A. Bartl, H. Fraas and W. Majerotto, Nucl. Phys. B278 (1986)1; G. Moortgat-Pick 
and H. Fraas, Eur. Phys. J. C 25 (2002) 189 | hep-ph/0204"333| . 

[8] M. Peskin, in Proc. of the Third Workshop on Physics and Experiments with Linear 
Colliders, ed. A. Miyamoto et al, (World Scientific, Singapore, 1996) p. 248; Int. J. 
Mod. Phys. C13 (1998) 2299 |hep-ph/9803279"l ; S. Thomas, Int. J. Mod. Phys. C13 
(1998) 2307 jhep-ph/9803"420| ; J.L. Feng and M.E. Peskin, Phys. Rev. D 64 (2001) 
115002 hep-ph/0105"T0"0"| ; A. Freitas, D.J. Miller and P.M. Zerwas, Eur. Phys. J. C21 
(2001) 361 |hep-ph/0106198| ; A. Datta, A. Djouadi and M. Muhlleitner, Eur. Phys. 
J. C25 (2002) 539 |hep-ph/0204354 ; C. Blochinger, H. Fraas, G. Moortgat-Pick and 



W. Porod, Eur. Phys. J. C24 (2002) 297 hep-ph/0201282|; J.A. Aguilar-Saavedra 

El 



and A.M. Teixeira, hep-ph/0307001 A. Freitas, A. von Manteuffel and P.M. Zerwas, 
|hep-ph/0310182 



[9] J.F. Gunion and H.E. Haber, Phys. Rev. D 37 (1988) 2515. 

[10] M.M. Nojiri and Y. Yamada, Phys. Rev. D 60 (1999) 015006 lhep-ph/990220l| |; A. 
Bartl, W. Majerotto and W. Porod, Phys. Lett. B465 (1999) 187 |hep- ph/9907377| ; 
A. Djouadi, Y. Mambrini and M. Muhlleitner, Eur. Phys. J. C20 (2001) 563 
|hep-ph/0104115| . 



13 



[11] H. Komatsu and J. Kubo, Phys. Lett. B157 (1985) 90; H.E. Haber and D. Wyler, 
Nucl. Phys. B323 (1989) 267; S. Ambrosanio and B. Mele, Phys. Rev. D 53 
(1996) 2541 |hep-ph/9508237] |; H. Baer and T. Krupovnickas, JHEP 0209 (2002) 38 
|hep-ph/0208277| . 

[12] M.S. Berger, Phys. Rev. D 41 (1990) 225; Y. Okada, M. Yamaguchi and T. Yanagida, 
Prog. Theor. Phys. 85 (1991) 1; Phys. Lett. B262 (1991) 54; J. Ellis, G. Ridolfi and F. 
Zwirner, Phys. Lett. B257 (1991) 83; H.E. Haber and R. Hempfhng, Phys. Rev. Lett. 
66 (1991) 1815. 

[13] H.E. Haber and Y. Nir, Phys. Lett. B306 (1993) 327; H.E. Haber, in Physics From 
the Planck Scale to the Electroweak Scale, Proc. of the US-Polish Workshop, Warsaw, 
Poland, September 21-24, 1994, edited by P. Nath, T. Taylor and S. Pokorski (World 
Scientific, Singapore, 1995) pp 49-63; J.F. Gunion and H.E. Haber, Phys. Rev. D 67 
(2003) 075019 |hep -ph/02070To"l . 

[14] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66 (2002) 01001. 



14 



